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ABSTRACT 

The planning for an underway at sea replenishment operation is 
treated as a special case of the job scheduling problem. Given the 
service times of the ships to be replenished, general expressions are 
obtained for total time to complete the operation and total waiting time 
of the ships involved. Methods to minimize these expressions are 
developed for several specific cases. Although the largest size of 
operation considered in this paper is that of three supply ships and 
three ships to be replenished the methods developed are believed to be 


extendable to larger operations. 
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I. INTRODUCTION TO THE PROBLEM. 

One of the techniques most heavily relied upon by the United States 
Navy in the projection of sea power is that of replenishment at sea. 

It is this technique which fs cut the strings of dependence on land 
bases which have been tied to combatant ships since steam replaced sail. 
Because of the assumption of global responsibilities by our nation, and 
sive vulnerability and political aspects of large fixed land bases, it is 
desirable to make naval task forces as independent of such fixed bases 
as possible. 

Both from a strategic and a tactical standpoint, replenishment at 
sea is one of the critical operations required of a naval force. While 
performing this operation, combatant ships are wedded to a relatively 
slow speed and a steady course. For the ships in any replenishment 
formation, there can be no immediate high-speed manuvers. The combat- 
ants are diverted from their normal operations to a situation which 

Maye s them more vulnerable to every form of attack than is desirable 
ome any length of time. In general, for the duration of the replen- 
ishment, the flexibility and responsiveness of the force is greatly 
RBeraded’ 

Eoaear cere at sea involves two important facets; planning for 
the operation, in a tactical sense, and the actual conduct of the 
operation. Considerable effort has been and is being expended to im- 
prove the latter. For example, individual ships strive to improve sea- 
manship and loading/unloading plans so as to decrease the time involved 
in replenishment. New types of supply vessels such as the USS SACRAMENTO 
(AOE-1) have been built in the past several years, combining the supply 
functions of two vessels in one ship. More efficient, computer assisted 


storage and unloading schemes have been instituted on some older and 


mewer vessels. The use of helicopters to deliver certain items, called 
vertical replenishment, has been developed, and with the continuing 
development of higher capacity helicopters such as the UH-46, one can 
envision the day when only the non-nuclear ships in a task force will be 
forced to conduct alongside replenishment for fuel (6). Unfortunately, 
problems such as the present limited strike-down capability of combatants 
are still to be resolved, and we must live for a while longer with the 
system we have. 

However, no matter how efficient the performance of individual units, 
the performance of the total operation is strongly dependent on the 
effectiveness of prior planning of the operation. The planning phase 
involves detailed knowledge of the capabilities of the units involved, 
and the interactions involved in the actual conduct of the operation. 
Certainly, one of the primary concerns of the planning process for 
current replenishment operations is the answering of the question "what 
ship goes where, and when?" It is a scheduling problem, and as currently 
done, is primarily based on experience and common sense. 

Analytical studies of the interaction and scheduling aspects of 
planning for an underway replenishment are few. McCullough (4) has 
studied the interaction aspect by considering the operation as a queueing 
process consisting of M service centers in series. The system serves N 
units, each of which goes through all stages in succession. In order to 
achieve a steady state solution, the units repeat the process continuously. 
He then applies the procedure developed by Koenigsberg for cyclic queues 
to investigate the effect of variation in service rate upon various 
system characteristics such as mean number of units waiting at a stage and 


time required for a unit to complete M stages. 


The most basic objection to this approach is that it is not a 
representation of what really happens in a replenishment at sea. First, 
the cyclic nature of this model is a great simplification of the more 
complex cycles involved in the operation. Second, the operation is 
finite, and therefore no steady state exists, forcing one to consider 
transient conditions. Last, there is the assumption by the author 
of exponential distributions on the service times to make the model more 
amenable to analysis. McCullough attempted to resolve some of these 
objections by developing a simulation program based on his model which 


more closely approximated reality. 


Il. OBJECTIVE 

An alternative approach which seems worthy of consideration comes 
from the scheduling area of the problem. This is the tack which we 
shall pursue in this paper. The objective shall be, in very general 
terms, to arrive at an analytical method for the determination of the 
"best'' schedule in given circumstances. More specifically, we will 
analyze the replenishment operation considered as a job scheduling 
problem where service times are assumed known. For the general case 
of m supply vessels and n combatants we will attempt to develop methods 
for determining the optimal solution to the replenishment scheduling 
problem for two criteria which shall be defined at a later point. 

It is appropriate at this point to mention possible bounds on the 
size of the problem in actual operations. To a fairly close approxi- 
mation, the largest size situation which might arise is on the order of 
m=6, n=20. An operation of this order of magnitude might be expected to 
occur in Sixth Fleet operations, but, as a matter of actual fact, in the 
more critical area of operations, the South China Sea, they rarely 
approach this size. In this area, because of various operational and 


tactical considerations, the size range should be closer tol2n2Z5, 


and 12 m= 4. 
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IIIT. REPLENISHMENT AS A JOB SCHEDULING PROBLEM 

If we let each of the m supply ships represent a machine, and 
each of the n combatants to be serviced represent a job, then our 
problem would be to find the optimal schedule for the n jobs on the m 
machines. The nature of underway replenishment introduces a unique 
twist to the job scheduling problem; that is the fact that several jobs 
may be started at once by different machines. A literature search re- 
vealed a considerable body of work done in the area of job scheduling! 
but none which was applicable to the problem with this twist. All of 
the work was concerned with the case in which all jobs must start at the 
first machine. Such a constraint would seldom, if ever be associated 
with replenishment at sea. Unfortunately, the removal of this constraint 
means that we must seek a general solution to a more complex combinator- 
ial problem than the "one line" case previously studied in the litera- 
ture. 

Elaborating on the notion of replenishment as a job scheduling 
problem, let us assume that we know the service time of each combatant 
at each supply ship, and that transit times between supply ships are 
identical for all combatants and negligible in comparison with service 
times. Further, let us assume that for this travel between supply ships, 
there exists a specified flow, say, from left to right for example. This 


means that the flow of jobs through the system can be represented as 





TSG, for example, Dudek and Teuton (1), Eastman, Even and Isaacs (2), 
Lomnicki (3), Palmer (5), and Smith and Teuton (7). 
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indicated in the figure below for the case m=3: 


D B b B'sie 


(\ COMBATANT 


Figure l 


For the moment we will assume that we have not specified a particular 
machine (supply ship) ordering. 

One may object that by assuming a specified flow pattern, we 
have deleted many possible solutions to the problem. While this is 
true, it is important to realize that the physical conditions under which 
a replenishment at sea operates preclude a large number of such possible 
flow patterns for reasons involving safety or good seamanship. It is a 
violation of the above principle to maneuver in the manner which would 
be required for certain flow patterns in the presence of the replenish- 
ment formation, which is a privileged formation. However, by assuming 
no specification of supply ship order, we do allow a larger number of 
patterns to be considered than appears feasible at first glance. 

The only feasible case not allowed by our formulation of the 


system is that which we will call "jockeying", described below for the 


Q0 Z Z 


Figure 2 


case m=3: 
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Suppose that, as depicted in Figure 2, the combatant on the left takes 
much longer than the right hand combination to complete its service; 
that is, the two combatants at the right may find themselves waiting 

for the left-hand combatant to finish its service. Then, the case may 
arise in which one would want to interchange the two rather than sending 
one or both to wait behind the left-most station. Note that the same 
situation may apply to any of the combinations of two stations in the 
system. For larger systems, the same type of situation may occur, but 


will be more complex. 
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IV. OPTIMALITY CRITERIA. 

There are two optimality criteria which appear appropriate for our 
problem. The first is the minimum total time to complete the entire 
replenishment operation, or "Min T" criterion. There is some schedule 
or set of schedules which will yield this Min T. This particular 
criterion may be the desirable path to follow when the operation is to 
be performed under a threat of some sort; air, surface, or sub-surface. 

The second criterion is less obvious. One might wish to minimize 
the total waiting time of all ships concerned, both combatant and 
supply. By definition, and in consonance with reality, a ship is not 
considered to be waiting once she has completed all her required services. 
This criterion might be applied where one is concerned more about the 
efficiency and smoothness of the operation than about total time to 
complete. We will call this the "Min M" criterion. This criterion 
will give total replenishment times which are greater than or equal to 
the value of Min T. 

It is important to note that careful consideration of the tactical 
situation will be necessary in making a choice between the criterion. 
For example, a situation may arise under threat for which it will be 
more desirable to use the Min M criterion to disperse the combatant 
force more rapidly, them to use the Min T criterion, which may keep a 


large number of combatants near the supply ships throughout the operation. 
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V. TOTAL WAITING AND COMPLETION TIMES 
A. Matrix of Service Times 
Let us assume there are m supply ships and n combatants par- 
ticipating in the replenishment operation. Let the supply ships be 
identified by the numbers one through m and the combatants by the 
numbers one through n. By letting Sij represent the known service time 
of combatant i at supply ship j, the nm services may be expressed as a 
matrix Ss » For example, in the case n = m = 3, 
ci = Ueeeris 
P= |so1 $22 893 
531 532 533 
The ith row of sy is the vector of service times for combatant i 
and the j column of S is the service vector of supply ship j. As 
defined this matrix of service times is dependent upon the indexing 
procedure used to number the combatant and supply ships. Since there 
are nim! ways in which the n combatants and m supply ships may be numbered 
there are necessarily n!m! unique matrices which could be constructed 
from the nm service times. The relationship between the n!m! matrices 
that could be obtained is that given any one we may generate the re- 
maining by elementary row and column transformations of the given matrix. 
In other words, there are n!m! ways in which the rows and columns of S 
may be arranged and each of these corresponds to one way of numbering 
the n combatants and m supply ships. The dependence of ss upon a 
particular numbering system can be removed by relaxing the usual 
definition for equality between matrices. For the purpose of this paper 
two service time matrices are equal provided they are row and/or column 


equivalent. 
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B. The Solution Matrix 

The flow or movement of combatants from one supply ship to the 
next has been assumed to be ordered and so a solution is complete if we 
merely specify the starting order. This specification will take the form 
of assigning stations to the m supply ships and starting positions to 
the n combatants. For the purpose of analysis no generality is lost by 
visualizing the m supply stations as being abreast and numbering the 
stations from left to right with the numbers one through m. When n = m 
the m combatant starting positions are taken to be those stations which 
initially are identical to the m supply stations. 

The question of combatant starting positions is not so easily 
handled when n is not equal to m. Such a case may be dealt with by 
creating either false supply ships or false combatants depending on 
whether n is greater than or less than m. A false ship of either type 
will have for its service vector a vector of zeroes. The use of false 
ships reduces all cases to the form m= n. To make this notion clear 


the case n = 3, m = 2 can be visualized as depicted in Figure 3. 


B B é 6) 


Y - SUPPLY 
STATION 


() _ COMBATANT 
STATION . 


Figure 3 
In this case one of the three supply stations is actually filled 
with a false supply ship which means that the combatant assigned to that 


station will be initially waiting for service. A combatant station 
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having once been assigned moves with the particular combatant assigned 
to that station. In the example above the positions after completion of 


the first services would be as given in Figure 4. 


@) @ fa) és) 


Figure 4 


Returning now to the matrix S , the elements S;;, may be given 


: 
a more specific meaning by the use of superscripts which represent their 
positions within the matrix and denote the assignment of stations. The 
element sty is thus the service time of combatant i at supply ship j 
and is also combatant station K at supply station L. With superscripts 
attached to its elements, each one of the (m!) 4 ways that Sg may be 
expressed becomes a solution to the scheduling problem. The main 
diagonal of a superscripted S matrix is the first set of services to 
be performed and this specification constitutes a solution. To determine 
an optimum solution, a rule must be devised that will indicate how to 
arrange the rows and columns of S so that when it is viewed as a 
solution matrix either total system completion time or total waiting 
time is a minimum. 

In devising such a rule it is convenient to associate a matrix 


A with the position of the elements SAB That is a 


— eb . 
j Sq: We will 


KL 
work henceforth with the position matrix A but it is important to 
understand that when we speak of minimizing of maximizing a function 


of the ai4'8: we are actually referring to the rearrangement of the rows 


and columns of g so as to achieve such a minimum or maximum. 
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C. Duplicate Solutions 


The matrix a! is unfortunately repetitious in solutions which 
can be obtaimed by the rearrangement of its rows and columns. For 
example if the rows and columns of the matrix {5 2 /7]| are rearranged 

3 3 4 
2G 


to either of the following forms: 4 1 2 36 4 eS no new 


information is obtained. The three matrices are identical solutions 
under our previous definitions. In the om case there are (mtye 
orderings of the rows and columns of S , but (m!)? (1-1/m) of these 
solutions are duplicates leaving only m!(m-1)! different solutions. To 
avoid such duplication we will always assign the supply ship having 
_ 2 5;; to supply station number 1. Such a convention does not 
further restrict our solution. 
D. Total Waiting Time 

The following notation will be used in the development of 

expressions for the waiting and completion times. If the operation 


starts at time zero, then: 


T = total system time, 
Wy = the time that the combatant assigned to station i completes 
its final service, 
t,; = the time that the supply ship assigned to station i completes 
its final service, 
in Fay, = the total of all service times for the combatant 


assigned to station i, 
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a.i = S ar; = the total of all service times for the supply ship 
c 


assigned to station j, 


= 
tl 


ij waiting time of the combatant assigned to station i for 


service at the supply ship assigned to station j. 


Wij waiting time of the supply ship assigned to station j for 
the combatant assigned to station i, 

Wi.= Swi; total waiting time for the combatant assigned to 
station i, 

W; = EM; = total waiting time for the supply ship assigned to 


station j, 


w= am = total waiting time of the combatants, 
L 

W = 2 W; = total waiting time of the supply ships. 
J 


Under the ordered flow which has been specified,several relations 


are immediately obvious: 


tE+1 ifn 
- le i=zan 


(2) (a) Wy = Ei saet 


Tee 


(b) Wy 


Also apparent is the fact that W = w, since: 


vi 
= Saar BG e Sty Seren allem Leo als Ba. 
L c l & C22 é é é 


> W 
The general expressions for Wij and Wij are: 
Lasit; is A es | “Wij-l ,0} 4 ? j 
(= 
J : 
0 Ll=J 


le 


ie em ao [5-4 eT era peek re 0 | 
1j 0 


where the subscript 0 = n and the subscript n+1 = 1. 


Now w and W may be written as, 


* 


mn nn ae nok n n 
SF; - 2(% Wjitj ap Sess] + ¢ ats Wiitj a ry nt 
| e 


W = 
U=t gst gin-trs Ce giNacr 
mon n-t n 

ie S FW; 5 - Abe ‘a rr r 5 iia + ON jj-nti 
ts| gr dine G41 cr) gil a= N-bei 


If we let M = W + w then, 


n-\ ne 
4 VAST + Wii+g) + SOs nes* en 


ttt y:' a: N-ctl 


and using equations (3) and (4), 


nes aie 
ca 2 (MAX [a jar ig jt pening? 2 jin g-1- jie j-10) + MAX Lagiyga + 
ea) aya) 


Wtj-1 Fylde 0]) 


[Me 


a P (vax Ays1jenti js jee 25 j-wti-1-"jj—mti- 10] ‘i 


a iam Or | 


MAX (85 joni -1TW 5 jembi-L = Qjat gonel = jet gened y al) ¢ 


a 


Finally, since MAX (A-B,0) + MAX (B-A,0O) = IA-BI, 


n-) ne 


TS y, S |= jrlitj Wjtlitj ?*jitj-17 Wyits- -1| i" S|: 4 j+1j-mei 


cst] yet yrn-t+y 





TW 54d j eh j j mt i-1 55 ont 


To minimize the waiting time within the system we have only to 
minimize M which is always twice the total waiting time. 
This fact is a direct result of the preceding work but can be 


justifiably questioned on the grounds of its interpretation. The use 
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of false ships necessarily causes their waiting times to be contained 

in M and these waiting times are, in a practical sense, of no concern 

to us. We are interested instead in the total waiting time of the ships 
that actually exist, the real total waiting time. Since we must either 
create false supply ships or false combatant ships, but never both, 
there is a simple relationship between total waiting time and real 

total waiting time. The quantity M/2, in terms of real waiting times, 


may be summarized as: 


The real total waiting time if false supply ships were 
of the combatants created, 

M/2 = }] The real total waiting time if false combatants were 
of the supply ships created, 
The real total waiting time if no false ships were 
of the combatants or supply created 
ships 


Notice that if a transit time, identical for all combatants, is 
added to each aij these constants would cancel in equation (5). We can 
thus ignore transit times provided they are assumed to be identical 
for all combatants. 

E. Total Completion Time 

It is apparent from equation (1) that the total system time T 
is given by, 

T = MAX olewaligs eee) = MAX (ty, t5,0tq.....t=)- 
Or from equation (2), (a) and (b), 


T = MAX (wy +a,.,Wotay....--Wi ta.) = MAX (Wytag .---- Wita gn) - 
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(6) T = MAX Sit, = MAX OW tal 


t 
FW 4+421 va Wjata.2 
P| L 


ar itn, 2 Winta n 


Having developed general expressions for M and T we will now turn 


our attention to their minimization. 


He 


VI. OPTIMAL SOLUTIONS 
A. Solutions for General m and n 

We have not as yet been successful in obtaining general solutions 
for the required minimizations of M and T. Our efforts to date have 
been directed toward ,solving specifie cases, i.e., given values of n, 
with the hope of eventually providing an inductive argument for the 
general case. Although it is unfortunate that no general solution has 
yet been found eicwmessite some consolation in the fact that the cases of 
small m and n are the ones that most often occur in the real world. 
Since the relevant cases are those of small m and n, optimal solutions 
for these few cases may be all that will ever be required. 

B. Solutions for Specific Cases 

The remainder of this report is concerned with the beginning 
of what we eventually hope will be a compilation of opt imal solutions 
for the relevant cases. In particular, solutions are presented for ° 
MIN M and MIN T in the cases n = 2 and n = 3. These solutions are in 
the immediately following ‘sections. . Section 5 presents a method- 
ology for obtaining MIN M by bounds. Although the procedure outlined 
is that for the case n = 3, we believe that the principle value of the 
method will come from its extension to larger cases. 

1. MIN M in the Case n = 2 

To minimize M we attempt to find some function of the aij's 

such that M in terms of this function will iridicate how the rows and 
columns of S should be arranged. In order to clarify this technique 
we will explore at some length the simple case n = 2. For this case 


equation (5) becomes: 


tas ylaggeerany iar” raza) ere ee wei 


igi 


First we note the obvious fact that the row and column sums of 
are arranged according to the way the rows and columns of S are 
arranged. This fact suggests that if M is expressed in terms of row 
and column sums and if this expression indicates how to arrange the 
row and column sums of S so as to minimize M, then we will also know 
how to arrange the rows and columns of pS) so as to minimize M. In 
terms of row and column sums, 

M = 2 lat - ay.| : 
Since 
gey = MAX=(S.),5h2)5 
in keeping with section C, it is obvious that to minimize M, 
ag. = MAX (S]_5S9_). 
We are told how to arrange the row and column sums of S and, therefore, 
how to arrange the rows and columns of S : 

A less obvious solution procedure is also worthy of consideration. 

We define 


5 Awe = Qe: - ae 


= 5. =n ij j ij 


i ple ig 
and construct the following matrix: 
7 Aiz> Aya Ata» Ayo 
Aoi, Agi Aga, Ago 
In these terms, 
M = 2 |Ag1 - Agil. 
which again suggests a method of arranging the rows of S : 
Consider the following example: 
g 4 3 = 552 Ri 
E S = 
2 1 4,1" 4258 


Since M is only a function of the pair of numbers that appear as the 


first element in the second row (the 21 position) of the 5 matrix, 
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it is apparent that to minimize M we must interchange the first and 
second row of Sg which corresponds identically with interchanging the 
first and second row of S : 

The number of other functions of the ai, s which could be used to 
solve for the minimum of M becomes larger as n increases. The problem 
is to find one for each case, or eventually, one for all cases which 
will lead to the optimum arrangement of the rows and columns of Se 

2. MIN M for the Case n = 3 
To obtain an expression for M in the case n = 3, equation 
(5) is expanded using the relation MAX (A-B,0) = MAX (A,B)-B. This 
gives: 
(7) Ms }211-422] + [222-433] + [233-411] 
+|a73+MAX (ap 9,433) ~a12-MAX (a1 1,299) 
+ | a31+MAX (a11 , 433) -a23-MAX(a22,233)] 
+| 41 +MAX (a) 1,499) -a37 -MAX (aj],233)| 
Equation (7) may be reduced using the relation )A-B] = 2MAX (A,B) -A-B and 


the definitions of Aij and A; to: 


j 

(8) M/2 = MAX(A,3,A13,A21,A32) + MAX(Ap1,A91 4325413) 

+ MAX (A39,A39 4732491) - (Ay 3tAg, +439) - 

This equation shows that M/2 depends only on the elements that appear 
in the 13, 21 and 32 positions of the Cl epri cs By the rearrangement of 
the rows and columns of a are twelve ways that these positions 
can be filled. It will be noted, however, that for any three given 
elements there are two ways in which they can fill the three positions. 
This observation permits us to divide the search for an optimum solution 
into two parts. We first identify the three elements from g which give 


rise to the optimum solution and then determine how these three elements 


are to be arranged to fill the 13, 21 and 32 positions of the matrix A. 
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To see how the second part is to be done, assume first that the 


= 
three elements that minimize M/2 are given from the rows of 5 SE 
these be: 
A13 =—8,5 A01 = 4,7 A320 = 6,6 


From equation (8) M/2 = 16+ 7 - 18 = 5 for this arrangement of 


the given elements, It is possible, however, by a rearrangement of 


=— 
= 


the rows and columns of S to obtain a second solution using these 
elements. That is: 

A13 = 6,6 Ag, = 4,7 R39 = 8,5 
For this arrangement M/2 = 16 + 6 - 18 = 4 which is the optimum solution. 
To determine which set of elements gives rise to the optimum solution 
reguires an inspection of each of the six choices available. For 
convenience the six choices may be expressed as the rows and columns 


of a matrix: 


= >in PEs aay 
S en 22 "BT ye 
ees = S21 | 


Since Aj 3+A9;+A39= Ai 344014439 it may be verified from eqeueeon (8) 
that the minimum value for M/2 for each row or column of S is 
determined by the following set of rules. 

1) If the MAX [Ais-A21 9A39,A13:491 »A30/ = Ai 


then, 
M/2 = A,, + MAX |Atct »Aam> Ag j MEN cy »Fgn) -Ax1 ~Aym 


k#i#n 
L# j#m 


2) Tithe MAX | Ats-A21-A32-A1 3.891 -Ag2]= Aiy 


then, r 
M/2 = Ay; + MAX |Aicd > Amm>Aa 5 sMIN Cte Avan) “Aye “Bien 


k#i#n 


1#j7m 
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A detailed example should serve to make the method clear. 


bet. Then : 
4 6 3 a Sa.V Te 7 gf ea 
@ = 7 ie S =/|7,8 12,8 9,4 
4 2 1 5,11 5,8 “Gx 
and; 
= oF ll 9,4 3,6 
S = {12,8 Pr 10,6 
6,8 7,4 7,8 


From the first column of g »M/2 = 12 + MAX |9, 6, 8,MIN(11,8)| 
- 9-6 = 6. 


By similar calculations the following results are obtained: 


> 9,11 9,4 5,8| M/2 = 8 
S = 112.8 3,11 10,6] M/2 =9 
6,8 7,4 7,8| M/2=4 

M/2=6 M/2=10 M/2=6 


Since the third row of Ss gives the minimum value for M/2, the elements 
which minimize M/2 are {6,8 7,4 7,83 . From a is apparent that 
by interchanging the second and third columns we will have: 

Ay, = 7,8 Ago = 6,8 and Ay3 = 7,4 


From equation (8) this arrangement gives: 


5 = 8+8+8- 20-4 
By appropriately interchanging the rows and columns of S we may also 
obtain: 

Aoi = 7,8 A329 = 7,4 and Ay3 = 6,8 


which gives M/2 = 4. There are thus two optimum solutions to the given 


problem. The first is obtained by interchanging the second and third 
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column of Ss and the second by interchanging the first and third rows 


of S . The optimal solutions are: 
4 3 6 Lima a 
S, = fer 5 Z and s = |/7 2 2) 
4 1 2 4 6 3 


For each solution the waiting time of the combatants is 4 and the 
waiting time of the supply ships is 4 since no false ships were created. 
3. MIN T for the Case n = 2 


For the case n = 2, equation (6) gives: 


T = MAX ay: 
Bo - 
aey 
Since 
MIN(MAX [ a,. |)= MAX a - 
a9 or 
| Gey 


any arrangement of the rows and columns of S is an optimum solution 
for T in the case n= 2, 
4. MIN T for the Case n = 3 
From equation (6), T for the case n = 3 becomes, after 
some reduction, 
T = MAX | MAX [293 + MAX(a99,433), ay + MAK (211>299)] + 443 
MAX [ a3] + MAX (aq1,233), a3 + MAX (a22,433)] + 221 


MAX | azo + MAX (2115499) » agi) lias (211,433) ] +. Bo 
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which may be expanded to: 


T = MAX ay. 


Tee 729 + 493 


491 + 493 + 433 





491 + 437 F 433 | 


Because in this equation MAX a,- | cannot be altered for any given 


oe 


ic 
ae] 
problem, every rearrangement of S will give the same value for this 


expression. This fact suggests that an alternative approach to seeking 


a solution having MIN T is to seek a solution giving a value of Z £O 


where, 
(9) 4A = MAX | aq] + a12 + 432 - MAX/a,. ; 
a faa] + aq0 a9. 
442 + 413 + 499 a3. 
a13 + a22 + 423 a-1| 


491 + 493 + 433 


a9, 7 437 = 32 


Any arrangement of S which yields a value of 4 €Q is an optimal 


solution for the minimization of T and therefore we need not seek that 


zo 


particular solution giving the MIN Z value. There could be perhaps 

a number of solutions for which ASO. This may be interpreted to 
mean that within the minimization of T there may exist some latitude 
for the adjustment of M. In equation (9) it may again be observed that 
the inclusion of transit times for combatants would result in their 
cancellation if all transit times are equal. An alternative form of 


equation (9) is given by equation (10). 


(10) A + MAX[ a). aj3 + a2 
= MAX 4192 + MAX 

ap « 411 + 439 

a3° 443 + 499 
293 + MAX 

aey a9] + 433 

ali 92 

a9) "486 


For an arrangement of S to be an optimum solution for MIN T the 


right hand side of equation (10) must be less than or equal to MAX ao 


The right hand side of (10) shows that the elements which appear in 
the 12, 23 and 31 positions of the Co matrix are the ones, subject to 
additive terms, that concern us. The term to be added to one of these 
position elements results from a matching of the elements in its row 


oi aa 
with those in its column. For the term aj + MAX ( ), Figure 


= Tilia 
5 shows the matching which pertains. 


30 


a) Fel 68 
“\S* 


= 


Figure 5 


For a + mat rix S there are twelve ways that the 12, 23 and 31 
positions may be filled. It may be shown, however, that T is only 
dependent on those numbers which fill these positions and not on the 
order in which they are filled. This means that there are at most six 
different T solutions (3!). For a given arrangement of the matrix Ss ; 
the six different choices to fill the three positions may be expressed 
as the rows and columns of a matrix A‘. The superscripts on the elements 


of A' denote only matrix position. 


11 12 pis 
“oT a39 13 
‘ 
= 21 22 23 
A ie ai Tal a59 
31 32 33 


ST? 993 431 
The values of aa for the six different choices of the elements to 
fill the 12, 23 and 31 positions of ert may be obtained by the use of 
the following equations. 
For the it? row of A’: 
ee ij Lk 
(11) = = MAX | atJimax ( a“){ - MAXfa,. : 
j k LF 
ao. 
ay. 


| eo 
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For the jth column of A‘; 
. ij a bik - ' 
OO ee ia ae MAX fa). 
i LL K# 
Ao. 
A3- 
a og, 


An example will show the procedure outlined above to be quite 


simple to apply. Consider the following matrix: 


Si thse 5 
Ce loon Cie 4 
ae 


2 5 § 

i] 
[oe 3 8 6 
30 


Since the MAX | a,.|= 17, equation (11) for the first row of A' gives 


NY 
it 


8+4 3+3 34+3 = 
max | 24max(€*9), Dass te? ieee ih | a 


MAX [ 15, 18,17] -17=1 


With equations (11) and (12) in mind the remaining cases may be 


examined without resorting to written form. Since MAX ay: is the 


a2 


same for each solution it is only necessary to scan a possible solution 
for combinations which disqualify it from giving a minimum value for T. 
For example, the second row of A' is disqualified by the combination 
8+5+7 = 20. By this procedure the only solution not disqualified is 
the first column of A' which gives % = 0. The elements of the first 
rwo of A' correspond to the elements in the 21, 33 and 12 positions of 
the matrix g : 

There are two ways to rearrange S so that these elements appear 
in the 21, 32 and 13 positions. The two solutions are: 


8 5 3 r/ 2 5 


= 
it 
NO 
> 
Ov 
NO 
tl 
NO 
Ov 
> 


7 3) 5 8 3 5 

Each of these solutions gives a value of T = 17 which is the 
minimum possible. 

There do exist cases for which MIN ~ > 0. These become obvious 
when it turns out that all possible solutions have been disqualified. 
In such a case it is only necessary to repeat the solution processes 
eliminating a particular solution only after a better solution is found. 


An example is: 


4 1 6 
ees ART oe 
5 6 1 
for which 
3 6 6 
‘ 
A = 1 4 Z 
1 4 5 


It is easily verified that MIN # = 1 and that there are four minimum 


solutions. 
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5. MIN M by Bounds for the Case n = 3 


There are twelve feasible solutions to our problem in the 
case n = 3. If upper and lower bounds can be placed on the value of M 
for these solutions then perhaps the set of solutions we must consider 
in detail may be a smaller set than the original twelve. We begin by 
determining the expressions for the upper and lower bounds on the value 
of M for a particular solution. 
From equations (1) and (2), 
Wy-ata.4-1 fl 
Wrtan i=1 
which can be written as: 


Ja.;.3-az-] #1 
[wi-Wi-1] i. 


ja.,-a,- | i=l 

The latter relation leads to a lower bound for M. Since there 
cannot exist a difference between waiting times unless waiting time 
exists, we can state that 

M 2 S [ aprai-| + [a-n-ay-| 


(= 9 
Using the notation: 


Range (A,B,C) R(A,B,C) the difference between the largest 


and smallest of A,B,C, 


Upper Range (A,B,C) =UR(A,B,C) 


the difference between the largest 
and second largest of A,B,C, 


and equation (7) we get the following upper bound for M when n = 3: 
M < 2UR(a11 5499433) + 2R(a,1,899,893) + 2R(a,9,893,€831)- 
The use of the upper and lower bounds on M to simplify the problem 


of finding MIN M can be shown most directly by considering a specific 


example. 
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Let 


4 E) 3 
le 
4 il! 2 

Then aj. is one of the set (14,10,7); a 15 and a., is one of the 


set (9,7). All possible pairings of the elements of these sets give 
the lower bounds for all the values of M. That is, we obtain the lower 
bounds 2, 6, 10, 12, 14 and 17. Starting with the pairing which gives 
the smallest member of this group, we calculate the upper bounds on M. 


For the lower bound M 8 2 there are two cases: 


ey < oe: 
S=- [x 2 5 ema = 7 eee 
4 2 @ y 3. = 


The upper bounds for 5 1 and Ses are 12 and 14 respectively. 

Since the upper bound on S 1 1s less than or equal to the lower 
bounds 12, 14 and 17, we may restrict our attention to the six distinct 
cases which give rise to lower bounds 2, 6 and 10. To decide between 
the cases remaining required a repeated application of equation (7). 

An alternate lower bound for M is: 

M2 2R (ay 1 +79 5433) +2R (Ay 7 »2935831) -2MIN] UR (ay 1 99,453) ,UR(217,8935831)] 
This inequality is easier to apply than equation (7) and may prove helpful 
in further reducing the set of cases that must be solved exactly by 
equation (7). In the example above it gives q 1 v7 8 which is the 
exact value for S 1’ 

, in this example the use of bounds on M reduced the problem by one 
hal €. Due to our lack of computational experience we cannot say that 
this amount of reduction will always occur whenever n = 3. The deter- 


mination of bounds for the values of M may somewhat simplify the problem 
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at hand. It is expected, however, that in the majority of cases such 
a reduction will still leave the problem unmanageable if we must cal- 
culate an exact solution for each of the remaining cases. Perhaps the 
most powerful use of this method is in its application as a test of 
optimality. For, if we can demonstrate one solution which has its 
value of M equal to the smallest Lower bound we know that solution is 


an optimal solution. 
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VII. DISCUSSION. 

We have achieved our objective only in a limited sense, in that 
we have developed solutions for an optimum schedule only in the cases 
m =n = 2 and m = n = 3. However, even these limited results are of 
value, for two reasons. First, as previously indicated, small system 
sizes of approximately this magnitude are the rule in at least one 
critical area of fleet operations, Southeast Asia. Second, the work 
we have done may point the way for further research which might lead 
to a general solution for any size operation. 

Larger problems will include even more complex combinatorial 
problems than the cases we have studied. An obvious approach towards 
solving large problems is to develop a computer program. Such an 
approach has been found practical even for the "one line" classical 
job scheduling problem when m 7 3 (1), (2), (3), (5), (7). However, 
we have avoided involvement with computer solution procedures because 
the majority of personnel involved with the replenishment scheduling 
function have neither a multipurpose computer nor a computer specialist 
available. Moreover, in trying to develop analytical solution techniques 
we have gained valuable insights into the replenishment operation which 
would have been lost in the inner workings of a computer had we taken 
this approach. 

One of these insights, in the form of a rather general statement 
concerning the effects of relatively large operation size on the choice 
of schedule, follows. As the operation becomes larger in the sense 
that the number of combatants involved begins to be much greater than 
the number of supply vessels involved, the planner has more latitude 


in his choices for optimum sequences. That is, he will be increasingly 


2/ 


less likely that any of the supply ships are not busy. The extent to 
which this effect occurs is a definite area for further research. 

There are several aspects of the actual underway replenishment 
operation which have not been included in the discussion, other than — 
perhaps in passing up to this point. The first of these is the technique 
of vertical replenishment. Since this technique is one which may well 
be a revolutionary factor in replenishment at sea in the future (and 
having already proven itself a useful adjunct to current operations), 
it should be mentioned. First, one should consider the question of the 
method of inclusion of vertical replenishment in the formulation of the 
problem as presented. Suppose that a particular combatant is capable 
of receiving all of the supplies it requires from a particular supply 
vessel by vertical replenishment. If this is the case, we propose that 
the vector of service times of the combatant in question include a zero 
as the service time of the combatant at the supply ship in question 
because the service can be completed while the combatant is either 
waiting for or completing another service. As a second proposal, we 
suggest that vertical replenishment be used as a method to make service 
times more even, that is, to decrease the spread of service times of 
the various combatants in order to make the system cycle more smoothly. 
Both of these suggestions are profitable areas for further study. 

It is a common practice in underway replenishment to assign one 
supply vessel to service two combatants at the same time. This is the 
second of the aspects of the operation which we have not taken into 
account. One possible method of handling this situation is to split 
each supply ship into two "effective" servicing units, one to represent 
the port side, and one the starboard side of the supply ship. Thus 


one might split a three by six problem into two separate three by three 
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cases, which can be handled by our methods. The port side of all supply 
vessels could represent one group, and the starboard side the other. 
Needless to say, there are many choices and criterion for the assign- 
ment of combatants to one group or the other, but, given an assignment, 
an optimum for the group is obtainable. Note that the spectre of sub- 
optimization rears its head in the consideration of this suggestion in 
that the assignment to a specific group must be done in an optimal 
manner in addition to optimal scheduling within the group, in order to 
insure an overall optimum. This area is likewise in need of consider- 
able further study. 

There exists a need for lifeguard vessels which we have not covered 
in our analyses to this point. As a safety precaution, some vessel 
must normally act as this lifeguard vessel for a replenishment station. 
That means that at any given time there should be m lifeguards on station, 
one for each of the m supply vessels. The usual procedure is to have 
some vessel of a smaller, more manuverable type, such as a destroyer, 
fill this station as it waits to be serviced. However, there is the 
possibility that there will be no such vessel available to perform the 
function. There are several methods by which one might handle this con- 
sideration if it appears that it will become a serious one. A group of 
m vessels of the required type may be detailed to act as permanent life- 
guards for each of the m supply ships, thus breaking the problem into 
two parts. The first would be the n - m combatant ships proceeding 
through the system first, and the second, them lifeguard ships. If 
there are enough helicopters available, they might be used as lifeguards. 
Lastly, if time is of sufficient criticality, one might ignore any 
problems of this sort which arise on the basis that the time considera- 


tion will outweigh the loss in safety. It should become obvious after 
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a schedule is established and checked just where and when such problems 
will arise. Decisions should be made at this point as to the most 
desirable remedy available. 

A- discussion of some of the more questionable assumptions we have 
made is now in order. First, and foremost, we have assumed that service 
times are known quantities. This assumption is obviously contrary to 
the true state of nature. The time involved has an unknown probability 
distribution varying generally with the amount of goods required, the 
type of goods required, with sea state and weather, and possibly many 
other factors. However, we feel that the determination of this distri- 
bution is possible, and is a necessary and important part of the develop- 
ment of optimum schedules in the real world. 

We envision the day when tables, or graphs, providing the re- 
quired service time information will be readily available to the planning 
agencies. If it can be shown that mean service times are sufficient, 


then perhaps graphs such as Figure 6 can be set up. 


Na 
MEAN oe 
TIME 
TO ° 
SERVICE Ss 


GALLONS REQUICED 


O FUEL OIilL~ DDX CLASS AT AOX CLASS 
(TRANSFER ZATE DEPENDENT ON TEMPE@ATURE) 


Figure 6 
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Currently, combatants inform tactical commanders and the replenishment 
group commander of their requirements for resupply by message, so that 
the information required by the abcissa is known. We envision a graph 
or curve of this sort for every class or type of combatant at each 

type of replenishment vessel. For the newer multipurpose supply ships 
such as the AOE and AFS classes, the information would have to be of 
such a nature as to take into account the diversity of supplies carried, 
as would the information for every type supply vessel except for oilers. 
The data necessary to build these graphs or curves exists, but consider- 
able effort would be required to sort out, collect, and analyze it. 

The assumption was made that the time spent in transit between 
supply vessels by combatants is both negligible in comparison to service 
times involved, and identical for all ships. We have shown that if they 
are identical then transit times cancel out and can be ignored in the 
actual application of our model. Therefore, the critical assumption 
is that the transit time is constant. While this is untrue, the vari- 
ation in these times is generally small enough in comparison with the 
service times, even for different type combatants, that the error intro- 
duced by the assumption will be slight. 

It has been stated that the use of computers to solve the problem 
of seater of replenishment at sea is not presently feasible due 
only to the lack of availability of either computers, or computer 
trained personnel, or both. As a proposal to alleviate this situation, 
we suggest the following innovation. The provision of a communications 
link which would allow the transmission of the required data, such as 
expected service times, from the tactical planning agency to a shore 
facility with a readily available computer and computer staff, should be 


investigated. Both the communications link and the computer service 
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should be of an on line nature. An arrangement of this sort could 
provide planning agencies with the necessary answers, and would allow 
the optimal solution of larger operations than is presently possible 


using our technique. 
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VIII. SUMMARY. 

In this paper we have investigated the optimal scheduling of 
underway replenishment operations from the starting point of consider- 
ing the problem as a job scheduling problem. We have developed general 
expressions for total waiting time of both combatants and supply vessels, 
and for the total completion time.:' Proceeding from this we have 
developed solution procedures for the optimal schedule under the two 
criteria Min M and Min T for the cases m= n= 2 andm=ne= 3. Al- 
though we have considered the cases individually, the possibility that 
there is an inductive connection between the first (smaller) cases and 
larger cases is not ruled out. 

The following areas related to the analyses in this paper present 
some of the most potentially profitable possibilities for further study: 

1. Extension of the applicability of analytical solution to 
a general solution. 

2. The extent of increase in latitude for schedule choice for 
large sized operations. 

3. The effect of the introduction of the proposed vertical 
replenishment modifications upon solutions and the operation 
as a whole. 

4, An investigation of the split supply ship concept. 

5. Investigation of the relationship between Min M and Min T 
solutions. 

6. A data collection and analysis effort in order to determine 
service time distributions. 

7. An investigation of the feasibility of the computer data 


link proposal. 
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